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FINAL  REPORT-  AOARD  064092 
"Target  Tracking  and  Interception  by  Aggressive  Honeybees" 


SUMMARY  OF  AIMS 

This  study  has  investigated  how  aggressive  bees  detect  and  track  their  targets,  to  uncover  how 
natural  visual  systems  have  evolved  to  perform  the  highly  exquisite  and  demanding  task  of 
detecting  a  moving  object  in  the  environment,  and  then  tracking  it  and  intercepting  it  quickly  and 
unerringly. 

SUMMARY  OF  RESEARCH  CONDUCTED 

We  have 

(i)  Designed  and  constructed  an  experimental  chamber  for  studying  and  filming  the 
behavior  of  aggressive  bees  (Fig.  1). 


Figure  1.  External  view  of  experimental  chamber,  showing  projector  and  computer  for  display 

of  targets  on  a  screen  inside  the  chamber. 
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ii)  Designed  and  developed  a  number  of  methods  for  initiation  of  aggression  in  honeybee 
colonies.  An  example  of  one  of  the  methods  is  shown  in  Fig.  2  and  is  described  later 
below. 


Figure  2.  View  of  specially  designed  “Aggravator”  to  create  aggression  in  a  honeybee  colony 
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Figure  3.  View  of  aggressive  bees  attacking  a  beekeeper’s  glove 


FINAL  REPORT-  AOARD  064092  -  M.V.  Srinivasan 


5 


(iii)  Set  up  and  calibrated  a  high-speed  stereo  camera  system  for  recording  the  flight 
trajectories  of  aggressive  bees. 


Figure  4.  View  of  inside  of  experimental  chamber,  as  captured  by  one  of  the  high-speed 
stereo  cameras  positioned  at  one  comer  of  the  chamber.  The  second  stereo  camera  (not 
visible)  is  in  the  ceiling.  The  hive  entrance  is  to  the  right. 

(iv)  Developed  computer  algorithms  for  digitizing,  tracking  and  reconstracting  the  body 
positions  and  orientations  of  individual  aggressive  bees  pursuing  targets  in  3  dimensions 
(see  Fig.  21  below). 

(v)  Compared  aggression  towards  stationary  3-dimensional  and  2-dimensional  targets. 

(vi)  Investigated  the  distribution  of  aggression  towards  multiple  stationary  targets. 

(vii)  Designed  and  tested  mechanically-based  as  well  as  video-display-based  equipment 
for  producing  moving  targets; 

(viii)  Set  up  a  computer-controlled  visual  display  to  investigate  how  aggression  depends 
upon  the  spectral  (colour)  properties  of  the  target 


(ix)  Acquired  substantial  quantitative  data  on  3-D  trajectories  of  aggressive  bees  pursuing 
moving  targets 
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(x)  Developed  a  preliminary  mathematical  model  to  characterize  the  dynamics  of  target 
pursuit  behavior  in  aggressive  honeybees 

(xi)  Investigated  a  novel  ‘streamlining’  response  displayed  by  tethered,  flying  bees, 
which  is  likely  to  be  relevant  to  achieving  high  flight  speeds  and  accelerations  during 
aggressive  pursuit. 


KEY  FINDINGS 

The  key  research  findings  are  summarized  below: 


1 .  Methods  of  inducing  aggression 

Two  methods  were  found  to  be  the  most  effective  in  eliciting  aggression  in  bee  colonies: 

(a)  Alarm  pheromone  method  Alarm  pheromone  was  obtained  by  removing  the  sting  and 
poison  sack  from  the  abdomens  of  a  few  bees  and  placing  them  on  a  piece  of  filter  paper,  which 
absorbed  the  extruded  pheromone.  The  filter  paper  was  then  placed  in  front  of  the  hive  entrance. 
After  about  a  minute,  aggressive  bees  emerged  from  the  hive  in  considerable  numbers.  A  movie 
illustrating  this  procedure  is  included  in  the  accompanying  Powerpoint  presentation. 

(b)  Mechanical  aggravator  method  A  specially  designed  device  consisting  of  an  array  of 
paint  brushes  mounted  on  a  sliding  rod  was  moved  back  and  forth  across  the  hive  entrance  (Fig. 
2).  The  moving  brushes  flicked  bees  emerging  at  the  entrance,  inducing  aggression  in  them.  A 
movie  illustrating  this  procedure  is  included  in  the  accompanying  Powerpoint  presentation. 

Depending  upon  the  particular  colony,  method  (a)  or  method  (b)  was  found  to  be  more  effective. 


2.  Experiments  to  investigate  target  preference 
(a)  3-D  versus  2-D  targets 

We  found  that  bees  in  a  state  of  aggression  will  track  and  intercept  three-dimensional  targets,  as 
well  as  two-dimensional  targets  projected  on  a  screen.  In  the  case  of  stationary  targets, 
aggressive  bees  prefer  to  attack  a  three-dimensional  target  (a  hemisphere),  as  opposed  to  a  two- 
dimensional  target  (a  flat  disc)  of  equal  surface  area. 

When  aggressive  bees  are  offered  a  choice  between  a  hemispherical  sphere  and  a  flat  disc  (of 
equal  diameter  or  equal  surface  area),  the  bees  display  a  greater  frequency  of  attacks  toward  the 
3-D  target  when  it  has  the  same  diameter  as  the  2-D  target,  but  a  lower  frequency  of  attacks 
when  it  has  the  same  surface  area.  This  was  established  by  video-filming  strikes  to  the  targets 
using  a  video  camera,  over  controlled  periods  of  time  following  initiation  of  aggression. 

These  findings  demonstrate  that  bees  are  capable  of  visually  distinguishing  between  thee- 
dimensional  and  two-dimensional  objects,  and  prefer  to  direct  their  aggression  at  the  three- 
dimensional  object  when  it  subtends  the  same  visual  angle  as  the  two-dimensional  object. 
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The  data  are  summarized  in  the  figures  below. 


Figure  5.  Total  number  of  strikes  to  a  2-D  target  and  a  3-D  target  placed  side  by  side.  The 
experiment  was  conducted  twice  (Series  1  and  Series  2)  and  repeated  each  time  by  swapping  the 
left-right  positions  of  the  two  targets,  to  average  out  any  side  biases.  The  data  represents  the 
pooled  results  of  experiments  with  3-D  and  2-D  targets  that  were  equal  in  (a)  diameter  or  (b) 
surface  area. 
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Figure  6.  Total  number  of  strikes  to  each  target  in  the  left  and  right  positions.  The  data  are 
pooled  across  replicates. 
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Figure  7.  Percentages  of  strikes  to  the  2-D  and  3-D  targets,  pooled  across  replicates  and  sides. 
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Figure  8.  Percentages  of  strikes  to  2-D  targets  that  are  of  the  same  diameter  (2Dd)  or  the 
same  surface  area  (2Dsa)  as  the  3-D  target.  Data  are  pooled  across  replicates  and  sides. 


(b)  Distribution  of  strikes  over  a  stationary  2-D  target 
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We  video-filmed  and  analyzed  the  spatial  distribution  of  strikes  over  a  stationary  2-D  target 
diameter  19  cm.  For  analysis,  the  target  was  subdivided  into  o  (a)  8  equal  45-deg  sectors,  and  (b) 
an  inner  circle  and  an  outer  annulus  of  equal  area,  as  shown  in  Fig.  9. 


Figure  9.  Illustration  of  subdivision  of  target  into  (a)  8  equal  45-deg  sectors,  and  (b)  an  inner 
circle  and  an  outer  annulus  of  equal  area.  In  (b)  the  radius  of  the  inner  circle  (r)  is  70.7%  of  the 
radius  of  the  target  (R),  which  ensures  that  the  inner  circle  and  the  annulus  are  of  the  same  area. 


The  number  of  strikes  to  each  of  these  regions  was  evaluated.  The  results  show  that  aggressive 
bees  prefer  to  strike  (a)  the  upper  half  of  the  target,  rather  than  the  lower  half  and  (b)  the 
boundary  of  the  target,  rather  than  the  interior.  The  strike  rate  is  highest  during  the  first  20 
seconds  of  aggression,  after  which  it  decays  approximately  exponentially  as  a  function  of  time. 
The  results  are  summarized  in  Figs.  10-13  below. 


Strikes  to  various  sectors  of  a  stationary  2D  target 


Sector  number 


Figure  10.  Strikes  to  each  of  8  sectors  of  a  stationary,  circular  2-D  target  of  diameter  19  cm.  For 
this  analysis,  the  target  is  subdivided  into  eight  equal  45-degree  sectors,  as  illustrated  in  Fig.  9. 
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Strikes  to  a  stationary  2-D  target  during  successive 
10-sec  intervals 


Figure  11.  Strikes  to  the  target  during  successive  10-second  intervals,  pooled  over  all  sectors. 


Number  of  strikes  in  the  inner  and  outer  target 
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Figure  12.  Comparison  of  strikes  to  the  inner  circle  and  the  annulus  of  the  2-D  target  (these 

areas  are  equal) 
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Number  of  strikes  in  the  inner  and  outer  target 

regions 


Inner  Circle  Outer  Circle 


Figure  13.  Total  numbers  of  strikes  to  the  inner  circle  and  outer  annulus  of  a  circular 
target,  both  regions  having  the  same  area.  Data  is  pooled  over  3  replicates. 


(d)  Attacks  to  multiple  targets 


To  examine  how  bees  direct  their  aggression  towards  multiple  targets,  we  presented  them 
with  a  constellation  of  9  dark,  circular  2-D  targets,  arranged  in  a  3x3  array  as  shown  in 
Fig.  14.  Strikes  to  the  targets  were  recorded  by  using  a  video  camera  to  fdm  the  display 
for  controlled  periods  of  time  following  initiation  of  aggression. 
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Figure  14.  Illustration  of  target  array  used  to  investigate  aggression  toward  multiple 
targets.  The  targets  were  each  9.4  cm  in  diameter,  and  were  separated  by  13  cm  (centre  to 
centre). 

The  results  (Fig.  15)  show  that  aggressive  bees  prefer  to  strike  the  peripherally  located 
targets  rather  than  centrally  located  ones.  The  central  target  elicits  the  lowest  strike  rate. 
Furthermore,  within  each  row  or  column  of  targets,  (or  across  the  two  diagonal  rows),  the 
middle  target  is  always  the  one  that  elicits  the  lowest  strike  rate.  This  may  reflect  a 
strategy  to  maximize  the  number  of  successful  strikes,  if  one  assumes  that  the  peripheral 
targets  can  hinder,  prevent  or  endanger  attacks  to  the  centrally-located  targets. 


Figure  15.  Percentage  of  strikes  (totalling  100%)  to  each  of  the  nine  targets  illustrated  in  Fig.  14. 

(c)  Spectral  properties  of  target  attractiveness 

We  find  that,  in  general,  aggressive  bees  prefer  to  strike  a  dark  target  on  a  light 
background.  We  have  commenced  investigating  the  spectral  properties  of  target 
attractiveness  by  comparing  the  strikes  elicited  by  white,  light  grey,  dark  grey  and  black 
targets  when  they  are  presented  against  a  green  background  or  a  blue  background,  as 
illustrated  in  Fig.  16.  All  of  the  targets  were  stationary.  The  positions  of  the  four  targets 
were  randomly  interchanged  from  trial  to  trial,  to  average  out  any  positional  biases  with 
regard  to  target  selection.  Strikes  to  the  targets  were  recorded  by  using  a  video  camera  to 
film  the  display  for  controlled  periods  of  time  following  initiation  of  aggression. 
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Figure  16.  Stimuli  for  investigating  spectral  properties  of  target  attractiveness. 

The  results  are  shown  in  Fig.  17.  In  general,  the  attractiveness  of  the  target  increases  as  it 
becomes  darker.  This  is  true  for  both  backgrounds,  except  for  the  one  case  of  the  black 
target  against  the  green  background.  A  black  target  against  a  blue  background  attracts 
the  largest  percentage  of  strikes.  Thus,  aerial  targets,  which  usually  appear  dark  when 
viewed  against  a  blue  sky,  are  likely  to  be  attacked  strongly,  as  are  terrestrial  targets 
when  they  are  approached  from  below  and  viewed  against  the  sky. 


Relative  preferences  for  targets 


Figure  17.  Relative  preferences  for  white,  light  grey,  dark  grey  and  black  targets 
presented  against  a  green  background  (green  curve)  or  a  blue  background  (blue  curve). 

(d)  Preference  for  stationary  versus  moving  targets 

Aggressive  bees  were  exposed  to  a  visual  display  projected  onto  the  inside  wall  of  the 
chamber  that  faced  the  hive  entrance.  The  display  presented  two  dark  targets  of  the  same 
diameter.  One  of  the  targets  was  stationary,  while  the  other  oscillated  from  side  to  side  at 
a  uniform  velocity  (target  position  versus  time  was  a  triangular  waveform).  Both  targets 
were  two-dimensional.  Strikes  to  the  targets  were  recorded  by  using  a  video  camera  to 
film  the  display  for  controlled  periods  of  time  following  initiation  of  aggression. 
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Figure  18.  View  of  video  camera  for  filming  and  recording  strikes  to  targets  and  background 
through  a  porthole  at  the  rear  of  the  chamber. 

The  results  showed  that,  during  the  first  minute  of  aggression,  bees  do  not  discriminate 
between  a  stationary  and  a  moving  target.  However,  they  show  a  clear  preference  for  the 
moving  target  during  the  second  and  third  minutes.  Thus,  bees  are  able  to  distinguish 
between  stationary  and  moving  targets,  whilst  they  themselves  are  in  motion.  The  data 
are  summarized  in  Figs.  19  and  20  below. 
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Moving  Static 


Figure  19.  Comparison  of  numbers  of  strikes  to  a  stationary  versus  a  moving  2-D  target. 


Moving/stationary  strike  ratio 


Figure  20.  Ratio  of  strikes  elicited  by  a  moving  versus  a  stationary  target,  measured  separately 
over  three  1 -minute  intervals  after  onset  of  aggression. 
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(2)  Tracking  and  interception  of  moving  targets 

We  have  filmed  a  large  number  of  flights  of  aggressive  bees  attacking  3-D  targets  that  (a) 
oscillate  along  a  straight  line  in  the  horizontal  plane  (b)  oscillate  along  a  straight  line  in  the 
vertical  plane  and  (c)  move  along  a  circular  trajectory  in  the  vertical  plane. 

An  example  of  a  pursuit  trajectory,  digitized  and  reconstructed  in  3-D,  is  shown  in  Fig. 

21.  An  animation  of  this  reconstruction  is  included  in  the  accompanying  Powerpoint 
presentation. 


Figure  21.  An  example  showing  the  reconstruction  of  the  3-D  trajectories  of  the  target 
(red)  and  a  pursuing  aggressive  bee  (orange).  The  red  dots  denoted  successive  positions 
of  the  centre  of  the  sphere,  and  orange  dots  and  lines  denote  successive  positions  and 
body  orientations  of  the  bee.  The  sphere  denotes  the  current  position  of  the  target. 

Preliminary  analysis  of  the  pursuit  trajectories  of  aggressive  bees  has  revealed  that: 

(a)  Pursuit  speed  increases  approximately  linearly  as  a  function  of  time,  i.e.  the 
acceleration  is  approximately  constant.  Two  examples  are  shown  in  Figs.  22  and 
23.  These  data  suggest  that  aggressive  pursuit  is  achieved  by  the  production  of  a 
constant,  high  thrust  and  that  the  acceleration  during  pursuit  is  limited  primarily 
by  the  inertia  of  the  bee,  and  not  by  aerodynamic  resistance. 
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Figure 


Bee  19 
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22. 


Flight  speed  versus  time  for  an  aggressive  bee  tracking  a  moving  target  (example  1). 


Bee  23 
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Figure  23.  Flight  speed  versus  time  for  an  aggressive  bee  tracking  a  moving  target 
(example  2). 

The  acceleration  continues  until  the  time  of  impact  with  the  target.  The  target  is  impacted 
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at  a  high  velocity  (about  2m/sec)  and  there  is  little  or  no  deceleration  prior  to  impact.  A 
movie  showing  the  impact  is  included  in  the  accompanying  Powerpoint  presentation. 

This  is  in  stark  contrast  to  the  landing  behavior  of  a  non-aggressive  foraging  bee,  which 
decelerates  smoothly  as  it  approaches  its  target  (a  food-reward-bearing  object).  Here,  the 
distance  to  the  target  (and  therefore  the  flight  speed)  decreases  exponentially  as  a 
function  of  time,  suggesting  that,  in  this  case,  the  landing  bee  adjusts  its  flight  speed  so  as 
to  hold  constant  the  rate  of  expansion  of  the  target  (Fig.  24).  Additionally,  foraging 
honeybees  display  a  hover  phase  just  prior  to  touchdown,  in  which  they  hover  at  a 
distance  of  about  15  mm  from  the  target.  This  hover  phase,  which  is  illustrated  in  a  movie 
in  the  accompanying  Powerpoint  presentation,  never  occurs  in  the  case  of  an  aggressive 
bee. 


S61 


Figure  24.  Distance  versus  time  for  a  foraging  (non-aggressive)  honeybee  approaching  a 
target  containing  a  food  reward.  The  red  circles  denote  mean  distance  as  a  function  of 
time,  the  green  curves  denote  ±  standard  deviation,  and  the  blue  curve  represents  a  fit  to 
an  exponential  function. 

(b)  the  turning  rate  of  the  bee  (yaw  rate)  is  controlled  by  the  angular  position  of  the 
target,  as  well  as  the  rate  of  change  of  the  angular  position.  This  is  illustrated  in 
Figs.  25  and  26. 
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Correlation  beeiotadion  rate  end  erroiarigle 


Figure  25.  Cross  correlation  between  bee  rotation  rate  (yaw  rate)  and  angular  position  of 
target  as  a  function  of  delay  time  (in  msec,  abscissa)  for  various  individual  bees. 

The  cross  correlation  between  the  bee  rotation  rate  and  the  angular  position  of  the  target 
is  a  maximum  when  the  target  angular  position  is  delayed  by  30  ms  relative  to  the  bee 
rotation  rate.  This  suggests  that  the  instantaneous  yaw  rate  of  the  bee  is  approximately 
proportional  to  the  angular  bearing  of  the  target  30  ms  earlier  in  time.  The  precise  value 
of  this  latency  varies  from  bee  to  bee. 
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Figure  26.  Cross  correlation  between  bee  rotation  rate  (yaw  rate)  and  rate  of  change  of 
angular  position  of  target  as  a  function  of  delay  (in  msec,  abscissa)  for  various  individual 
bees. 

The  cross  correlation  between  the  bee  rotation  rate  and  the  rate  of  change  angular 
position  of  the  target  is  a  maximum  when  the  rate  of  change  of  angular  position  is 
delayed  by  approximately  0  ms  relative  to  the  bee  rotation  rate.  This  suggests  that  the 
instantaneous  yaw  rate  of  the  bee  is  approximately  proportional  to  the  current  rate  of 
change  of  the  angular  bearing  of  the  target. 

On  the  basis  of  the  above  data,  we  have  formulated  a  preliminary  model  to  characterize 
the  dynamics  of  pursuit  behavior  of  aggressive  bees  tracking  a  moving  target.  The  yaw 
rate  co  of  the  pursuing  bee  is  given  by 

o)  =  a6{t-T)  + p^t)  (1) 

where  0  is  the  angular  bearing  of  the  target  in  the  bee’s  eye,  and  ¥  is  the  rate  of  change 
of  this  angular  bearing  (see  Fig.  27).  T  is  the  time  delay  associated  with  the  response  to 
the  target  bearing,  a  is  a  gain  parameter  that  relates  the  yaw  rate  to  the  angular  bearing  of 
the  target,  and  p  is  a  gain  parameter  that  relates  the  yawjj'ate  to  the  rate  of  change  of  the 
angular  bearing  of  the  target. 
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Target 


Figure  27.  Illustration  of  dynamical  model  for  characterizing  pursuit  behavior.  6  is  the 
angular  bearing  of  the  target  in  the  bee’s  eye,  and  ^  is  the  rate  of  change  of  this  angular 
bearing. 

We  have  carried  out  preliminary  simulations  of  this  model,  which  we  have  implemented 
in  Matlab.  One  example  of  the  results  is  shown  in  Fig.  28,  which  compares  the  observed 
trajectory  of  a  bee  (black  symbols)  pursuing  a  target,  with  the  trajectory  predicted  by  the 
model  (red  symbols),  using  the  following  parameters: 

(a)  Constant  acceleration  of  4000mm/sec 

(b)  a=20.0 

(c)  p=0.7 

(d)  T  =0.03  sec 

It  is  evident  that  the  model  captures  the  essential  elements  of  the  dynamics  of  the  pursuit 
behaviour.  An  animation  of  the  simulation  is  included  in  the  accompanying  Powerpoint 
presentation. 

Further  work  is  required  to  fine-tune  the  mode  and  to  evaluate  its  applicability  to  the 
large  number  of  pursuit  events  that  we  have  recorded. 
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Figure  28.  Illustration  of  one  example  of  the  results  from  a  model  simulating  the 
dynamics  of  pursuit  behavior.  The  circles  show  successive  positions  of  the  target  [blue 
(past)  red  (present)].  The  black  symbols  denote  the  reconstructed  trajectory  of  an 
aggressive  bee  pursuing  the  target,  and  the  red  symbols  denote  the  trajectory  predicted  by 
the  model. 

3.  Investigation  of  a  ‘streamlining’  response  in  flying  honeybees 

During  the  course  of  this  research  we  have  become  interested  in  examining  how  honeybees 
achieve  the  high  flight  speeds  and  accelerations  that  we  have  observed  when  they  are  in  the 
aggressive  state.  We  have  recently  discovered  a  novel  ‘streamlining’  response,  exhibited  by 
tethered,  flying  honeybees,  in  which  the  abdomen  is  held  in  a  raised  position  when  the  visual 
system  is  exposed  to  a  pattern  of  image  motion  (i.e.  a  flow  field)  that  is  characteristic  of  forward 
flight.  We  believe  that  this  visually  evoked  response,  which  can  be  elicited  without  exposing  the 
insect  to  any  airflow,  serves  to  reduce  the  aerodynamic  drag  that  would  otherwise  be  produced 
by  the  abdomen  during  real  flight.  The  streamlining  response  is  critically  dependent  on  the 
presence  of  appropriate  image  motion  in  the  frontal  as  well  as  the  rear  fields  of  view.  It  is 
relatively  insensitive  to  variations  of  stimulus  contrast,  exhibiting  values  close  to  the  maximum 
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for  contrasts  in  the  range  of  0.3  -  1 .0,  and  is  therefore  capable  of  operating  robustly  in  a  wide 
variety  of  natural  environments.  These  findings  also  underscore  the  importance  of  using 
panoramic  stimulation  for  the  study  of  visually  guided  flight  in  insects. 

The  experimental  procedures  and  key  results  are  summarized  below,  and  a  movie  illustrating  the 
‘streamlining’  response  is  included  in  the  accompanying  Powerpoint  presentation. 
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Figure  29 


The  abdominal  posture  of  tethered  honeybees  was  examined  as  they  flew  in 
a  virtual  reality  setup.  The  bees  were  tethered  to  the  horizontal  end  of  an  L- 
shaped  metallic  rod.  Flight  was  initiated  by  lowering  a  supporting  platform. 
{G]  Each  honeybee  was  placed  In  the  centre  of  an  arena  composed  of  four 
monitors.  A  video  camera  was  used  to  Him  a  side  view  of  the  bee  in  flight.  ^C) 
The  virtual  scene  consisted  of  an  infinite  tunnel  composed  of  two  walls  lined 
with  vertically  oriented  red  and  white  sinusoidal  gratings.  The  honeybee  was 
at  a  virtual  distance  of  X  units  from  each  of  the  tunnel  walls,  and  expenenced 
virtual  forward  motion  at  a  speed  of  Y  units  per  display  frame. 
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Figure  ^0 

Illustration  of  the  measurement  of  the  streamlining  response,  defined  as  the 
orientation  of  the  abdomen  relative  to  the  thorax.  The  upper  row  of  sketches 
illustrates  the  definition  of  the  response ,  and  the  lower  row  shows  two 
examples  of  how  this  is  used  to  measure  the  response. 
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Figure  31 

Comparison  of  streamlining  responses  Induced  by  forward  image  motion 
{FIM),  no  image  motion  (NiM)  and  backward  image  motion  (BiM).  The  grey 
staircase  displays  the  progressively  ascending  stimulus  speed  (see 
Methods).  Each  data  point  shows  mean  abdomen  angie  ±  standard  enor  of 
the  mean.  Note  the  comparatively  small  variability  in  the  responses  when 
bees  were  experiencing  forward  image  motion. 
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Figure  32 

Streamlining  responses  of  tethered  bees  in  experiments  using  various 
Donfiguralions  of  two  monitors  {coioured  symbols),  compared  with  the  four- 
monitor  conliguralion  (black  symbols).  All  flights  were  presented  wilh  the 
ascending  image  motion  speed  protocoi,  as  indicated  by  the  grey  staircase. 
Each  data  point  shows  mean  abdomen  angie  ±  standard  deviation.  Oniy  data 
from  bees  that  repeatedly  exhibited  continuous  flight  were  anaiysed. 
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Figufie  33 


Effect  of  stimulus  contrast  on  the  streamlining  response.  Variation  of  abdomen 
orientation  with  stimulus  speed  for  five  different  contrasts,  as  shown  in  the 
index.  Each  data  point  shows  the  mean  abdomen  angle  ±  standard  error  of  the 
mean. 
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A  manuscript  describing  this  study  is  in  preparation  (copy  attached). 

Although  AFOSR  funding  for  this  project  has  now  ceased,  we  will  continue  to  work  on  this 
project  with  the  aid  of  other  funding  sources.  We  will  continue  to  acknowledge  AFOSR  support 
in  research  manuscripts  currently  under  preparation  as  well  as  future  manuscripts  describing  the 
results  of  work  arising  from  that  initially  funded  by  AFOSR. 

We  wish  to  express  our  gratitude  to  the  AFOSR  for  their  initial  funding  in  making  this  study 
possible. 
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35  Summary 

36 

37  Animals  that  move  rapidly  through  the  air  can  save  considerable  energy  by  reducing  the  drag 

38  that  they  need  to  overcome  during  flight.  We  describe  a  novel  ‘streamlining’  response  in 

39  tethered,  flying  honeybees  in  which  the  abdomen  is  held  in  a  raised  position  when  the  visual 

40  system  is  exposed  to  a  pattern  of  image  motion  that  is  characteristic  of  forward  flight.  This 

41  visually  evoked  response,  which  can  be  elicited  without  exposing  the  insect  to  any  airflow, 

42  presumably  serves  to  reduce  the  aerodynamic  drag  that  would  otherwise  be  produced  by  the 

43  abdomen  during  real  flight.  The  response  is  critically  dependent  on  the  presence  of 

44  appropriate  image  motion  in  the  frontal  as  well  as  the  rear  fields  of  view.  Thus,  our  results 

45  also  underscore  the  importance  of  using  panoramic  stimulation  for  the  study  of  visually 

46  guided  flight  in  insects. 

47 

48 

49 

50 
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5 1  Introduction 

52  It  is  abundantly  clear  that  all  creatures  that  move  rapidly  —  in  the  air,  on  the  ground  or  under 

53  water  -  would  save  considerable  locomotive  energy  if  their  bodies  were  shaped,  or  adjusted, 

54  to  minimize  the  drag  that  is  produced  by  the  medium  through  which  they  move.  For  example, 

55  sharks,  tuna  and  mackerel  possess  fusiform  shapes  that  enable  high  swimming  speeds  to  be 

56  achieved  with  relatively  low  drag  (Schmidt-Nielsen,  1984;  Bailey,  1997). 

57 

58  While  it  is  of  interest  to  examine  whether  the  bodies  of  certain  animals  have  evolved  to 

59  reduce  drag,  it  is  also  pertinent  to  ask  whether  animals  can  actively  change,  or  adapt,  the 

60  shapes  of  their  bodies  to  cope  with  the  increased  drag  that  inevitably  accompanies  higher 

61  speeds  of  locomotion.  This  would  seem  to  be  important,  because  the  force  due  to  drag 

62  theoretically  increases  as  the  square  of  the  speed  at  which  the  body  moves  through  the 

63  medium  (Batchelor,  2000).  In  principle,  an  animal  could  achieve  ‘active’  streamlining  of  this 

64  kind  by  sensing  its  motion  through  the  medium  either  mechanically  through  the  use,  of,  say, 

65  mechanosensory  hairs  or  organs  to  sense  airspeed  during  flight  (Goodman,  2003),  or  through 

66  vision,  by  sensing  the  motion  of  the  image  of  the  environment  in  the  eye. 

67 

68  Do  animals  achieve  active  streamlining?  Here  we  ask  whether  honeybees  use  vision-based 

69  cues  to  streamline  their  bodies  at  high  flight  speeds.  We  find  that  bees  indeed  possess  a 

70  visually-driven  ‘streamlining’  response,  whereby  an  increase  in  the  image  velocity 

71  experienced  by  the  bee  causes  the  abdomen  to  be  raised  higher,  thus  partially  counteracting 

72  the  increased  drag. 

73 

74  As  an  insect  flies  through  the  environment,  its  eyes  experience  various  patterns  of  image 

75  motion  -  also  known  as  ‘optic  flow’  patterns  -  which  depend  upon  how  the  insect  moves  in 

76  relation  to  its  surroundings.  Research  over  the  past  fifty  years  has  shown  that  insect  flight  is 

77  controlled  by  optic  flow  in  a  number  of  different  ways.  For  example,  the  well-known 

78  ‘optomotor’  response  enables  a  flying  insect  to  stabilize  its  attitude  with  respect  to  the  world 

79  by  sensing  the  patterns  of  optic  flow  that  are  generated  by  its  rotations  about  the  yaw,  pitch 

80  and  roll  axes  (or  combinations  thereof)  and  delivering  appropriate  commands  to  the  flight 

81  motor  system  to  counteract  these  rotations  (Gotz,  1965;  Reichardt,  1969;  Srinivasan  and 

82  Bernard,  1977;  Collett  et  ah,  1993;  Krapp  and  Hengstenberg,  1996;  Egelhaaf  et  al.,  2004). 

83  Fruit  flies  (David,  1979,  1982;  Fry  et  al.,  2009)  and  honeybees  (Srinivasan  et  al.,  1996;  Baird 

84  et  al.,  2005)  regulate  the  speed  of  their  flight  by  holding  constant  the  magnitude  of  the  optic 
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85  flow  that  is  experienced  by  the  eyes.  Collisions  with  obstacles  are  avoided  by  detecting  and 

86  steering  away  from  regions  in  the  visual  field  that  are  associated  with  unusually  large  image 

87  speeds  (Srinivasan  et  ah,  1993;  Srinivasan  and  Zhang,  1997)  or  rapid  expansion  of  the  image 

88  (Tammero  and  Dickinson,  2002).  Collision-free  flight  through  narrow  passages  is  achieved 

89  by  balancing  the  magnitudes  of  the  optic  flow  that  are  experienced  by  the  two  eyes 

90  (Srinivasan  et  ah,  1991;  Srinivasan  et  ah,  1993;  Srinivasan  et  ah,  1996).  Smooth  landings  on 

91  a  horizontal  surface  are  achieved  by  holding  constant  the  rate  of  motion  of  the  image  of  the 

92  surface  as  the  surface  is  approached  (Srinivasan  et  ah,  2000). 

93 

94  Here  we  examine  whether,  and,  if  so,  how,  optic  flow  affects  body  posture  during  flight,  by 

95  filming  tethered  bees  flying  in  a  virtual-reality  arena  that  simulates  the  optic  flow  that  would 

96  be  generated  by  forward  flight  in  a  tunnel.  We  find  that  the  simulated  forward  flight  evokes  a 

97  ‘streamlining’  response  in  which  the  bee  raises  its  abdomen,  lifting  it  progressively  higher  as 

98  the  speed  of  the  image  is  increased.  This  streamlining  response  is  mediated  visually,  because 

99  it  is  displayed  by  tethered  bees  flying  in  still  air,  and  ceases  when  the  image  stops  moving. 

100  Presumably,  the  streamlining  response  serves  to  reduce  the  aerodynamic  drag  that  would 

101  otherwise  be  produced  by  the  abdomen  during  real  flight. 

102 

103 

104 

105  Materials  and  Methods 

106 

107  Bees:  Adult  forager  honeybees  were  used  in  all  of  the  experiments.  They  were  identified  as 

108  those  carrying  pollen  on  their  hind  legs.  The  bees  were  anesthetized  in  a  refrigerator  for  20- 

109  30  min,  after  which  they  were  taken  out  one  at  a  time  for  tethering. 

110 

111  Tethering:  While  the  bee  was  anaesthetized,  the  base  of  an  L-shaped  metal  rod  (see  Fig.  1) 

112  was  attached  to  the  thorax  by  a  small  globule  of  dental  glue,  which  was  cured  using  blue 

113  (440-480  nm)  light  (~40  s).  Adhesion  to  the  tether  was  facilitated  by  gently  shaving  the  hair 

114  on  the  thorax  using  a  scalpel  blade.  The  bee  was  then  housed  in  a  Styrofoam  box  in  which  the 

115  temperature  was  maintained  at  26-28°C.  A  beaker  of  water,  placed  inside  the  box,  provided 

116  an  environment  with  the  appropriate  humidity. 

117 
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118  The  experimental  arena:  Four  24”  LCD  monitors  (Dell  2408WFP)  were  arranged  on  a  table 

119  to  form  a  square  arena.  A  small  (~2  cm)  gap  was  provided  at  the  two  lateral  comers  of  the 

120  arena  to  allow  a  video  camera  to  film  the  bee  in  side  view  (Fig.  IB).  The  bee  was  tethered 

121  and  positioned  in  the  centre  of  the  virtual  environment  by  attaching  the  upper  end  of  the  L- 

122  shaped  rod  to  an  adjustable  metallic  arm  (not  shown).  This  arm  allowed  the  head  of  the  bee  to 

123  be  positioned  in  the  centre  of  the  arena,  at  a  height  midway  between  the  base  and  top  of  the 

124  monitors.  The  bee  faced  one  comer  of  the  arena.  The  table  was  covered  with  a  60  cm  x  90  cm 

125  sheet  of  black  cardboard.  All  experiments  were  performed  under  fluorescent  lighting. 

126 

127  Virtual  reality  environment:  Images  were  created  on  the  four  monitors  to  provide  a 

128  panoramic  virtual  environment  (Fig.  1).  The  LCD  monitors  were  driven  by  a  PC  (Intel  Quad 

129  Core  CPU  Q9300  2.5GHz,  3GB  RAM,  Windows  XP  2002,  DirectX  9.0)  with  two  dual  head 

130  NVIDIA  GeForce  8800  GT  video  cards.  The  monitors  had  the  following  settings:  resolution 

131  1920x1200  pixels,  refresh  rate  60Hz.  Since  the  monitors  used  LCD  displays,  there  was  no 

132  60Hz  flicker. 

133 

134  A  program  was  written  that  simulated  the  visual  effects  of  the  movement  of  an  observer 

135  within  a  custom-designed  virtual  world.  The  virtual  world  consisted  of  an  (apparently) 

136  infinite  tunnel  with  a  blue  ceiling,  a  dark  floor,  and  side  walls  lined  with  a  sinusoidal  grating 

137  composed  of  alternating  red  and  white  bars.  The  red  bars  would  have  appeared  as  (almost) 

138  black  bars  to  the  bee,  because  the  colour  red  falls  outside  the  bee’s  visible  spectmm  (Menzel 

139  and  Blakers,  1976;  Chittka  et  al.,  1993).  Red  was  used  rather  than  black,  to  facilitate  viewing 

140  and  filming  of  the  tethered  bee  against  the  background  of  the  visual  display. 

141 

142  [Insert  Fig.  1  here] 

143 

144  The  program  to  create  the  display  was  written  in  C++  and  used  the  DirectX  API  software  to 

145  create  four  ‘viewports’  of  the  virtual  world,  each  displayed  by  one  monitor,  with  the  virtual 

146  camera  position  geometrically  matched  to  the  centre  of  the  arena.  Each  monitor  displayed  one 

147  quadrant  of  the  virtual  world  -  forward  right,  forward  left,  rear  left,  and  rear  right.  A  software 

148  interface,  written  in  Python,  enabled  customization  of  the  virtual  world  and  control  of  the 

149  virtual  position  and  motion  of  the  observer  within  it.  For  a  bee  flying  along  a  virtual  tunnel, 

150  the  maximum  image  velocity  as  experienced  by  the  eyes  (the  maximum  optic  flow  rate,  OFR) 

151  occurs  in  the  lateral  viewing  direction,  i.e.  in  a  viewing  direction  at  90°  to  the  direction  of 
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flight.  For  a  simulated  linear  grating  speed  of  Y  units/ffame,  a  frame  rate  of  60/s  and  a 
simulated  flight  along  a  trajectory  that  is  X  units  from  either  wall,  this  maximum  OFR  (in  °/s) 
is  given  by 


OFR{°/s) 


60FA 


^180^ 
y  TT  J 


We  used  a  stimulus  protocol  that  simulated  flight  at  a  progressively  increasing  speed,  namely, 
Y  =  0.02,  0.03,  0.04,  0.05,  0.06,  0.07,  0.08,  0.09  and  0.1  units/frame.  The  width  of  the 
simulated  tunnel  (2X)  was  0.8  units.  The  height  of  the  simulated  tunnel  was  1.0  unit.  This 
staircase  of  simulated  flight  speeds  generated  the  following  OFRs:  172,  258,  344,  430,  516, 
602,  688,  773  and  859  °/s  respectively.  Each  epoch  of  stimulus  speed  was  8s  in  duration,  thus 
producing  a  stimulus  that  lasted  for  72  s  in  each  trial. 

Calibration  of  stimulus  contrast 

Sinusoidal  gratings  of  four  different  contrast  levels  were  used.  These  contrasts  (as  defined 
below)  were  nominally  set  to  1.0,  0.75.  0.5  and  0.4,  for  both  the  green  and  blue  channels  in 
the  24-bit  RGB  bitmaps  generated  using  Matlab.  The  red  channel  was  kept  constant  at 
maximum  value  over  the  entire  bitmap.  These  grating  bitmaps  were  tessellated  to  build  the 
long  virtual  tunnels  (described  above).  However,  subsequent  measurement  and  calibration  of 
the  various  gratings,  as  described  below,  revealed  that  the  actual  contrasts  of  these  gratings 
were  0.72,  0.54.  0.36  and  0.29  (or  about  72%  of  the  nominal  contrast).  It  is  the  actual  values 
that  are  used  to  report  the  results. 

The  contrasts  of  the  gratings  were  measured  as  follows.  First,  the  radiant  spectrum  of  the 
monitors  was  measured  using  an  Ocean  Optics  USB2000+  spectrometer  calibrated  for 
wavelength  and  intensity.  Spectra  were  measured  for  illumination  from  a  small  patch  within 
the  centre  of  the  white  bar  of  the  sinusoidal  grating,  as  well  as  within  the  centre  of  the  red 
bar.  The  relative  photon  flux  emitted  by  each  of  these  patches  was  obtained  by  integrating 
their  respective  spectra  over  a  spectral  window  of  300  nm  -  600  nm,  which  represents  the 
visible  spectrum  of  the  honeybee  (Menzel  and  Blakers,  1976;  Chittka  et  al.,  1993).  Denoting 
the  results  of  this  integration  for  the  white  and  the  red  bars  by  Imax  and  Imin,  respectively,  the 


6 


184  contrast  C  of  the  grating  was  calculated  as  C  =  .  This  quantity, 

185  known  as  the  Michelson  contrast,  is  the  commonly  used  measure  of  contrast  for  sinusoidal 

186  gratings  (Michelson,  1927).  It  can  assume  a  maximum  value  of  1.0  when  the  dark  bar  emits 

187  no  light  (Imin  =  0),  and  a  minimum  value  of  0.0  when  the  dark  and  bright  bars  have  the  same 

188  intensity  (Uax  =  Imin)- 

189 

190  In  all  experiments  but  one,  the  Michelson  contrast  of  the  grating  was  set  at  a  high  value 

191  (C=0.72  in  the  experiments  shown  in  Figs  3,  4  and  5).  In  one  experiment,  which  investigated 

192  the  effect  of  stimulus  contrast,  response  versus  speed  curves  were  obtained  for  a  series  of 

193  contrasts,  as  shown  in  figure  6. 

194 

195  Data  acquisition  and  analysis 

196 

197  A  Samsung  video  camera  (VP-HMX20C,  25  frames  per  second,  1920x1080  resolution)  was 

198  used  to  film  the  side  view  of  the  bee.  The  video  clips  were  analyzed  to  measure  the 

199  orientation  of  various  parts  of  the  body,  using  a  Matlab  digitization  program  written  in-house. 

200  The  program  computed  the  orientation  of  the  line  joining  any  two  manually  digitized  points 

201  in  the  image.  It  was  used  to  measure  the  orientation  of  two  lines.  One  line,  connecting  the 

202  head  to  the  thorax,  was  defined  to  be  the  orientation  of  the  thorax.  Another  line,  connecting 

203  the  front  and  rear  ends  of  the  abdomen,  was  defined  to  be  the  orientation  of  its  long  axis.  The 

204  program  calculated  the  difference  between  these  two  orientations  (Fig.  2).  This  difference, 

205  representing  the  orientation  of  the  abdomen  relative  to  the  thorax,  was  defined  as  the 

206  abdomen  angle,  or  the  response.  The  response  was  positive  or  negative,  according  to  whether 

207  the  abdomen  was  elevated  or  depressed  relative  to  the  thorax.  Responses  measured  in  this 

208  way  were  averaged  over  the  last  4  seconds  of  each  8-second  epoch  of  stimulus  speed,  thus 

209  giving  an  average  measurement  of  the  honeybee’s  steady-state  response  to  each  of  the  9 

210  stimulus  speeds.  Typically,  the  abdominal  response  reached  a  new  steady  state  within  l-3s 

211  following  each  new  stimulus  speed  (data  not  shown). 

212 

213  Statistical  analysis 

214 

215  All  statistical  analyses  were  performed  using  the  software  package  SPSS/PASW  Statistics 

216  18.0  (2009).  We  analyzed  the  abdominal  responses  from  two  independent  series  of 
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217  experiments  to  evaluate  the  effect  of  grating  contrast  and  monitor  configuration,  respectively. 

218  In  each  case,  we  first  performed  two-way  repeated  measures  ANOVA  to  investigate  the 

219  effect  of  each  independent  variable,  at  all  OFRs  of  the  stimulus  protocol.  If  a  significant 

220  effect  was  detected,  post-hoc  pairwise  comparisons  (t- tests,  with  Bonferroni  adjustment  for 

221  the  family- wise  Type  I  error  rate)  were  performed.  Note  that  for  the  mains  effect  of  interest  in 

222  each  case  (monitor  configuration  or  grating  contrast),  there  was  no  detectable  violation  of  the 

223  sphericity  assumption. 

224 

225  Comparison  of  responses  elicited  by  two  and  four  monitors:  Abdominal  responses  were 

226  compared  across  the  various  two-monitor  configurations,  and  across  two-  and  four-monitor 

227  configurations,  all  at  full  grating  contrast  (C=0.72).  Only  responses  from  honeybees  that 

228  experienced  all  the  different  configurations  were  used  for  this  particular  comparison  (n=7). 

229  Since  a  significant  effect  was  only  detected  when  the  4-monitor  data  was  included,  post-hoc 

230  tests  were  only  performed  in  the  latter  case.  These  statistical  results  are  shown  in  Tables  A1 

231  andA2. 

232 

233  Variation  of  response  with  stimulus  contrast:  Abdominal  angles  were  compared  across  the 

234  various  stimulus  contrast  levels  (0.29,  0.36,  0.54  and  0.72),  with  and  without  the  control 

235  condition  of  zero-contrast  (0.00),  all  using  the  four-monitor  configuration.  Only  responses 

236  from  honeybees  that  experienced  all  the  different  contrasts  (including  0.00)  were  used  for  this 

237  particular  comparison  (n=6).  Since  a  significant  effect  was  only  detected  when  the  zero- 

238  contrast  data  was  included,  post-hoc  tests  were  only  performed  in  that  case.  These  statistical 

239  results  are  shown  in  Tables  A3  and  A4. 

240 

241 

242  [Insert  Fig.  2  here] 

243 

244  Flight  initiation 

245  For  each  flight  trial,  a  tethered  honeybee  was  placed  in  the  centre  of  the  arena  and  a  platform 

246  (Fig.  1)  was  raised  from  under  the  bee  until  all  of  the  bee’s  tarsae  made  contact  with  it,  with 

247  the  legs  in  their  natural  resting  position.  As  soon  as  the  visual  stimulus  was  started,  the 

248  platform  was  retracted  from  under  the  bee.  At  the  end  of  the  visual  stimulus,  the  platform 

249  was  restored  to  its  pre-flight  position  so  that  the  bee’s  legs  were  again  resting  on  it,  thus 

250  preventing  further  flight.  The  platform  was  covered  with  wax  paper  (which  was  frequently 
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251  replaced)  to  ensure  a  firm,  but  smooth  support  that  could  not  be  gripped  by  the  bee’s  tarsal 

252  hooks,  thus  ensuring  trouble-free  release  of  the  feet  at  the  onset  of  flight.  After  each  flight  the 

253  bee  was  fed  10  pL  of  1  M  sucrose  solution  and  returned  to  the  humid  box. 

254 

255  Results 

256 

257  We  found  that  in  all  bees  tested,  the  action  of  removing  the  platform  from  under  the  bee’s 

258  legs  immediately  triggered  wing  beat  activity.  This  is  the  well-known  ‘tarsal  reflex’  (Fraenkel 

259  1932;  Chadwick  1953).  At  this  instant,  the  abdomen  would  also  drop  and  assume  a  hanging 

260  orientation.  Bees  would  then  occasionally  fly  for  a  few  seconds.  However,  combining  the 

261  platform  removal  with  the  onset  of  a  moving  visual  stimulus  enhanced  the  duration  of  flight 

262  considerably,  typically  at  least  two  minutes,  as  we  shall  describe  below.  Hence,  our  protocol 

263  was  to  place  the  tethered  honeybee  in  the  centre  of  the  arena,  rest  it  comfortably  on  the 

264  platform,  and  then  commence  the  stimulus  motion  as  soon  as  possible  after  platform 

265  withdrawal  (typically,  within  a  second). 

266 

267  Honeybees  only  exhibit  flight  behaviour  during  simulated  forward  motion 

268 

269  Withdrawal  of  the  platform  always  triggered  wing  beat  activity.  However,  this  activity 

270  typically  ceased  seconds  thereafter  if  there  was  no  image  motion  (NIM).  The  abdomen  would 

271  then  hang  almost  vertically  and  often  curl  inwards  toward  the  head  (Fig.  3).  Similar 

272  responses  were  observed  when  honeybees  were  exposed  to  backward  image  motion  (BIM), 

273  namely,  image  motion  that  simulated  backward  flight.  In  contrast,  even  on  the  very  first  test 

274  flight,  the  majority  of  honeybees  (76  out  of  101)  that  experienced  forward  image  motion 

275  (FIM)  -  image  motion  that  simulated  forward  flight  —  exhibited  continuous  flight  for  the 

276  entire  duration  of  the  visual  stimulus.  Presentation  of  FIM  in  the  two  front  monitors  and  BIM 

277  in  the  two  rear  monitors  -  or  vice  versa  -  also  evoked  no  flight  (unpublished  observations). 

278 

279  A  novel  ‘streamlining  ’  response  observed  during  simulated  forward  image  motion 

280 

281  Bees  presented  with  forward  image  motion  responded  not  only  with  wing  beat  activity,  but 

282  also  with  a  ‘streamlining’  response.  As  the  speed  of  the  image  motion  was  increased,  the 

283  position  of  the  abdomen  would  rise  progressively  and  then  remain  in  the  elevated  position  for 
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284  the  rest  of  the  trial  (see  Fig.  3,  and  movie  in  Supplementary  Material).  The  abdomen  reached 

285  its  highest  elevation  (approximately  -10°)  at  an  image  speed  of  430  °/s  and  remained  at  this 

286  elevation  for  all  higher  speeds,  up  to  the  maximum  speed  of  859  °/s.  Thus,  at  speeds  beyond 

287  430  °/s,  the  streamlining  response  was  maximal  and  approximately  independent  of  stimulus 

288  speed. 

289 

290  [Insert  Fig.  3  here] 

291 

292  Panoramic  motion  vision  in  honeybees 

293 

294  Bees  were  also  tested  in  the  same  arena  with  different  combinations  of  only  two  active 

295  monitors  (the  other  two  monitors  were  switched  off,  creating  a  black  display).  All  tests  were 

296  conducted  using  the  same  stimulus  protocol  of  ascending  image  speed.  The  bees  exhibited  a 

297  streamlining  response  even  with  the  two  monitor  setups.  However,  in  this  case  they  flew  less 

298  consistently,  and  for  shorter  periods.  Fig.  4  shows  the  relationship  between  the  response  and 

299  stimulus  speed  for  various  combinations  of  stimulation  with  two  monitors:  (i)  left  and  right 

300  front  fields,  (ii)  left  and  right  rear  fields,  (iii)  left  front  and  right  rear  fields,  and  (iv)  right 

301  front  and  left  rear  fields.  Statistical  analyses  were  performed  to  see  whether  the  abdominal 

302  responses  for  all  of  these  four  different  two-monitor  configurations  differed  significantly.  The 

303  two-way  repeated  measures  ANOVA  did  not  reveal  any  significant  mains  effect  of  monitor 

304  configuration,  F(3,  18)=0.44,/)=0.73  (see  Table  Al). 

305 

306  However,  it  is  clear  that  in  comparison  with  the  results  obtained  using  four  monitors  the 

307  responses  obtained  with  two  monitors  are  consistently  lower.  This  is  also  reflected  in  the 

308  statistical  results.  We  find  that  there  is  a  significant  mains  effect  of  monitor  configuration, 

309  F(4,  24)=19.62, /»<0.001  (see  Table  A2),  when  the  dataset  includes  responses  from  four  and 

310  all  two-monitor  configurations.  Post-hoc  pairwise  comparison  of  the  monitor  configurations 

311  at  each  OFR  showed  significant  differences  between  the  four-monitor  configuration  and  each 

312  two-monitor  configuration  ,  (p<0.005  with  Bonferroni  correction  in  all  cases,  see  Table  A2). 

313 

314  Furthermore,  while  the  curve  for  four  monitors  is  essentially  flat  for  speeds  greater  than  430 

315  °/s,  the  two-monitor  curves  exhibit  a  clear  peak  at  stimulus  speeds  between  602-773  °/s  (Fig. 

316  4).  With  two  monitors,  similar  response-versus-speed  curves  are  obtained  regardless  of  the 

317  configuration  of  the  monitors  -  all  four  configurations  yield  similar  results  (Fig.  4). 
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318 

319  [Insert  Fig.  4  here] 

320 

321  One  should  eonsider  the  possibility  that  the  intermediate  peak  observed  in  the  two-monitor 

322  tuning  curves  arises  as  a  result  of  fatigue  during  flight  through  the  prolonged  (72-second) 

323  stimulus  protocol,  causing  the  response  to  diminish  at  the  higher  speeds.  This  is  unlikely, 

324  however,  because  the  phenomenon  is  not  observed  with  four-monitor  stimulation. 

325  Nevertheless,  we  performed  control  experiments  with  two  monitors  in  which  the  stimulus 

326  speed  was  not  increased  progressively,  but  instead  held  constant  for  the  entire  period  of  72 

327  seconds.  Under  these  conditions,  we  found  that  the  response  was  nearly  constant  throughout 

328  the  period  of  stimulation.  This  is  shown  in  Fig.  5  for  two  different  stimulus  speeds.  These 

329  results  demonstrate  that  the  bell  shaped  speed  tuning  curves  obtained  with  two-monitor 

330  stimulation  are  not  a  consequence  of  response  fatigue  -  they  represent  a  variation  of  the 

331  sensitivity  of  the  response  to  the  speed  of  the  stimulus. 

332 

333  [Insert  Fig.  5  here] 

334 

335  Effect  of  stimulus  contrast 

336  How  sensitive  is  the  streamlining  response  to  variations  in  the  contrast  of  the  visual 

337  environment?  This  was  investigated  by  conducting  experiments  using  the  standard  stimulus 

338  protocol  (as  in  the  experiments  of  Figs  3  and  4)  for  different  Michelson  contrasts  (C)  of  the 

339  grating.  The  results  are  shown  in  Fig.  6.  The  profiles  of  the  curves  for  C=0.72,  C=0.54, 

340  C=0.36  and  C=0.29  are  very  similar  in  shape  and  height. 

341 

342  Statistical  analyses  were  performed  to  examine  whether  the  abdominal  responses  for  all  of 

343  these  four  grating  contrasts  differed  significantly.  These  analyses  did  not  detect  any 

344  significant  mains  effect  of  contrast,  F(3,  15)=1.47,  p=0.26  (see  Table  A3).  Consistent  with 

345  the  ANOVA  results,  the  corrected  post-hoc  pairwise  comparisons  between  the  non-zero 

346  contrast  conditions  did  not  show  any  significant  differences  (in  all  cases,  p>0J5).  On  the 

347  other  hand,  when  the  abdominal  responses  obtained  from  the  same  bees  for  a  grating  contrast 

348  of  0.0  were  included  for  comparison,  results  showed  that  there  was  a  significant  mains  effect 

349  of  contrast,  F(4,  20)=53.24,  /?<0.001.  Post-hoc  pairwise  comparison  of  contrast  conditions  at 

350  each  OFR  showed  significant  differences  between  each  zero  contrast  and  all  corresponding 

351  non-zero  contrast  conditions  (in  all  cases, /)<0.005  with  Bonferroni  correction,  see  Table  A4). 
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352 

353  This  means  that  the  streamlining  response  goes  from  zero  to  a  value  that  is  close  to  its 

354  maximum  value  when  the  contrast  of  the  grating  is  increased  from  0.0  to  0.29,  and  then 

355  remains  more  or  less  constant  for  contrast  values  ranging  from  0.29  to  0.72.  In  other  words, 

356  the  streamlining  response  is  relatively  robust  to  changes  in  contrast  in  the  mid- to  high 

357  contrast  range  of  0.29  -  0.72. 

358 

359  [Insert  Fig.  6  here] 

360 

361  Discussion 

362 

363  Our  results  have  revealed,  for  the  first  time,  the  existence  of  a  visually-driven  ‘streamlining 

364  response’  in  flying  insects.  This  response,  which  can  be  elicited  without  exposing  the  insect 

365  to  any  airflow,  presumably  serves  to  reduce  the  aerodynamic  drag  that  would  otherwise  be 

366  produced  by  the  abdomen  during  real  flight.  While  earlier  studies  with  flying  insects  have 

367  documented  elevation  of  the  abdomen  during  flight  (Moth:  Willmott  and  Ellington,  1997; 

368  Bumblebee:  Dudley  and  Ellington,  1990;  Drosophila:  David,  1978;  Honeybee:  Nachtigall  et 

369  ah,  1971;  Esch  et  ah,  1975),  they  observed  this  elevation  in  free,  un-tethered  flight.  Hence,  it 

370  is  difficult  to  ascertain  from  these  studies  whether  the  raising  of  the  abdomen  was  meditated 

371  by  the  visual  input,  by  mechanosensory  signals  driven  by  the  wind,  or  simply  by  passive 

372  mechanical  raising  of  the  abdomen  by  the  wind.  However,  Zanker  (1988)  reports  a  cursory 

373  observation  of  a  small  difference  (±2°)  in  abdominal  deflection  between  backward  image 

374  motion  and  forward  image  motion  in  Drosophila.  Our  study  demonstrates  comprehensively 

375  that  the  streamlining  response  can  be  elicited  by  a  visual  stimulus  alone,  namely,  a  moving 

376  pattern,  and  documents  its  properties. 

377 

378  The  streamlining  response  is  highly  selective  to  the  direction  of  motion  of  the  visual  pattern. 

379  Forward  image  motion  (motion  in  the  direction  that  would  be  experienced  during  forward 

380  flight)  creates  a  pronounced  elevation  of  the  abdomen.  On  the  other  hand,  a  stationary 

381  pattern,  or  a  pattern  that  moves  in  the  opposite  direction  (simulating  backward  flight)  does 

382  not  elicit  any  response  -  the  abdomen  remains  in  its  lowered  state  (Fig.  3). 

383 

384  Panoramic  stimulation  with  image  motion  seems  to  be  critical  for  eliciting  a  strong 

385  streamlining  response.  With  four  monitors,  the  response  is  significantly  and  substantially 
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386  stronger  than  that  measured  with  two  monitors,  at  all  stimulus  speeds  (Fig.  4).  The  four- 

387  monitor  stimulus  produces  a  streamlining  response  that  increases  steadily  with  stimulus  speed 

388  up  to  a  speed  of  about  400  °/s,  and  then  remains  constant  at  this  maximum  value  for  all 

389  higher  stimulus  speeds.  The  two-monitor  stimulus,  on  the  other  hand,  yields  a  bell-shaped 

390  curve  for  response  versus  stimulus  speed,  exhibiting  a  peak  at  stimulus  speeds  in  the  range 

391  602-773  °/s,  regardless  of  the  particular  configuration  of  the  two  stimulating  monitors  (Fig. 

392  4). 

393 

394  These  findings  reveal  that,  when  studying  visually  evoked  responses  in  flying  insects,  it  is 

395  important  to  ensure  that  the  stimulation  covers  all,  or  as  much  of  the  visual  field  of  the  insect 

396  as  possible.  Panoramic  visual  stimulation,  which  is  what  an  insect  would  experience  during 

397  actual  flight  in  a  real  environment,  reveals  response  properties  that  are  very  different  from 

398  those  produced  by  visual  stimulation  of  restricted  regions  of  the  visual  field.  Consequently, 

399  many  of  the  earlier  studies  of  visually-evoked  behaviour  in  insects  that  did  not  use 

400  panoramic,  or  nearly  full-field  stimulation,  would  have  to  be  reassessed. 

401 

402  What  causes  four-monitor  stimulation  to  be  more  effective  than  two-monitor  stimulation?  A 

403  simple  possibility  is  that  the  responses  of  movement-detecting  neurons  that  sample  various 

404  sub-regions  of  the  bee’s  panoramic  visual  environment  are  spatially  summed  to  generate  a 

405  command  signal  that  drives  the  streamlining  response.  Stimulation  of  the  visual  system  with 

406  four  monitors,  rather  than  two,  would  then  produce  a  stronger  streamlining  response.  This 

407  hypothesis  of  summation  is  also  consistent  with  our  observation  that  presentation  of  FIM  in 

408  the  two  front  monitors  and  BIM  in  the  two  rear  monitors  -  or  vice  versa  -  also  evokes  no 

409  flight  (data  not  shown),  presumably  because  the  oppositely-directed  stimuli  cancel  each  other 

410  out,  in  each  case. 

411 

412  If  such  spatial  summation  does  indeed  occur,  the  summation  must  not  be  linear  because  the 

413  response-speed  curve  that  is  obtained  with  four  monitors  caimot  be  predicted  accurately 

414  simply  by  scaling  the  two-monitor  curve  by  a  factor  of  two.  However,  the  four-monitor  curve 

415  can  be  explained  if  we  postulate  that  the  streamlining  response  is  driven  by  a  saturating 

416  summation  of  the  responses  of  the  various  movement-detecting  neurons.  The  possibility  is 

417  supported  by  the  observation  that  the  standard  deviations  of  the  responses  obtained  with  the 

418  four-  monitor  stimulus  are  consistently  smaller  than  those  obtained  with  two  monitors,  at 
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419  each  of  the  tested  speeds  except  for  the  lowest  speed  (172  °/s)  -  a  property  that  is 

420  characteristic  of  saturation  at  high  response  levels. 

421 

422  Can  the  saturation  of  the  streamlining  response  be  prevented  by  lowering  the  contrast  of  the 

423  visual  stimulus  to  reduce  the  strength  of  the  movement  signal?  Evidently  not.  The  data  in  Fig. 

424  6  show  that,  with  the  four-monitor  stimulation,  the  response  curves  obtained  for  grating 

425  contrasts  of  0.54,  0.36  and  0.29  are  remarkably  similar  to,  and  not  significantly  different  from 

426  that  obtained  for  a  grating  contrast  of  0.72.  This  must  mean  that  the  responses  of  the 

427  movement  detecting  neurons  that  drive  the  streamlining  response  are  relatively  insensitive  to 

428  variations  of  stimulus  contrast  in  the  range  0.29  -  0.72.  This  in  turn  implies  that  this 

429  invariance  to  contrast  must  come  about  through  an  invariance  of  each  movement-detecting 

430  neuron’s  response  to  the  stimulus,  for  contrasts  exceeding  about  0.3.  The  lamina  neurons, 

431  which  are  second-order  neurons  of  the  insect  visual  pathway,  exhibit  response  saturation  at 

432  contrasts  exceeding  0.4  (Laughlin,  1981).  If  these  lamina  neurons  provide  the  inputs  to  the 

433  movement-detecting  neurons,  their  properties  would  partially  account  for  the  postulated 

434  contrast  invariance.  Indeed,  individual  movement-detecting  neurons  in  the  insect  visual 

435  pathway  show  response  invariance  at  contrasts  exceeding  about  0.3  (e.g.  Dvorak  et  al.,  1980). 

436  Thus,  all  of  the  properties  of  the  streamlining  response  that  have  been  observed  so  far  can  be 

437  explained  by  postulating  phenomena  occurring  at  two  levels  of  the  visual  pathway:  (a) 

438  contrast  invariance  of  the  responses  of  the  movement-detecting  neurons  and  (b)  saturating 

439  summation  of  the  responses  of  these  movement-detecting  neurons. 

440 

441  The  finding  that  the  streamlining  response  already  attains  a  strength  close  to  its  maximum 

442  value  for  sinusoidal  gratings  with  a  Michelson  contrast  as  low  as  0.29  -  which  corresponds  to 

443  a  root-mean-square  (RMS)  contrast  of  0.20  -  is  very  compatible  with  the  fact  that  natural 

444  environments  generally  present  RMS  contrasts  in  the  range  of  0.2-0.3  (e.g.  Frazor  and 

445  Geisler,  2006).  Thus,  the  characteristics  of  the  visually  driven  streamlining  response  are 

446  tailored  to  the  properties  of  the  natural  environment. 

447 

448  The  study  by  Nacthigall  and  Hanauer-Thiesser  (1992),  measuring  the  aerodynamic  drag  that 

449  is  experienced  by  honeybees  for  various  pitches  of  the  abdomen,  demonstrates  clearly  that 

450  drag  is  reduced  when  the  abdomen  is  raised,  and  is  minimal  when  the  abdomen  is  oriented 

451  horizontally,  i.e.,  parallel  to  the  direction  of  flight,  when  it  presents  a  minimal  cross-section 

452  to  the  horizontal  airflow.  Thus,  the  changes  that  we  observe  in  abdominal  pitch  could  indeed 
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453  reflect  a  visually  driven  response  that  serves  to  reduce  drag  as  the  perceived  speed  through 

454  the  environment  is  increased.  Indeed,  freely  flying  bees  tend  to  hold  the  abdomen  high  when 

455  cruising,  and  low  when  hovering  or  flying  at  slow  speeds  (Esch  et  al.,  1975).  In  real  flight, 

456  however,  bees  should,  at  least  in  principle,  be  able  to  estimate  their  flight  speed  not  just 

457  through  visual  cues,  but  also  through  mechanosensory  cues.  For  example,  changes  in 

458  abdominal  pitch  could  be  mediated  by  a  reflex  that  is  driven  by  sensing  airspeed  as  signaled 

459  by  the  Johnston’s  organ  at  the  base  of  the  antennae,  and/or  other  wind-sensitive  hairs  on  the 

460  head  and  body  (Goodman,  2003).  Additionally,  the  abdomen  could  simply  be  pushed  up 

461  passively  at  higher  flight  speeds  by  the  increased  aerodynamic  pressure  acting  on  the 

462  underside  of  the  abdomen  and  the  legs  (Nachtigall  et  al.,  1971;  Nachtigall  and  Hanauer- 

463  Thiesser,  1992).  Notwithstanding  these  other  possibilities,  our  experiments,  conducted  with 

464  tethered  bees  in  still  air,  demonstrate  clearly  that  vision  plays  an  important  role  in  tuning  the 

465  posture  of  the  body  during  flight. 

466  Nachtigall  and  Hanauer-Thiesser  (1992)  have  shown  experimentally,  by  blowing  wind  at 

467  various  velocities  against  the  trunk  of  a  honeybee,  that  the  drag  force  increases  as 

468  (velocityf  ^  This  implies  that  when  a  bee  increases  its  flight  speed  from  a  value  that  is  close 

469  to  hover  (say,  5  cm/sec)  to  a  value  that  is  typical  of  cruise  (say,  10  m/sec)  the  drag  would 

470  increase  by  a  factor  of  about  2,800,  if  the  body  did  not  change  its  posture  or  shape.  The  same 

471  study  finds  that  decreasing  the  angle  of  attack  of  the  body  from  20°  to  0°  can  reduce  the  drag 

472  by  about  50%,  which  would  lead  to  a  significant  increase  of  flight  speed  and/or  a  significant 

473  saving  of  energy  over  a  long  flight.  Thus,  it  would  appear  that  actively  raising  the  abdomen 

474  indeed  produces  a  significant  and  beneficial  streamlining  effect. 

475  Finally,  the  variation  of  posture  with  flight  speed  may  represent  more  than  just  a  streamlining 

476  response-  it  may  reflect  the  fundamental  means  by  which  the  insect  adjusts  its  thrust  and  lift 

477  to  achieve  flight  at  various  speeds.  Horizontal  flight  at  higher  speeds  requires  the  wings  to 

478  produce  not  only  increased  thrust,  but  also  decreased  lift  -  because  the  abdomen  and  the  legs 

479  generate  greater  passive  aerodynamic  lift  at  higher  speeds  (Nachtigall  and  Hanauer-Thiesser, 

480  1992).  It  has  been  suggested  that,  like  the  fruitfly  Drosophila  (Gotz  and  Wandell,  1984; 

481  Zanker,  1988;  Vogel,  1996),  honeybees  are  unable  to  vary  the  direction  of  the  net  force 

482  vector  that  is  generated  by  their  wings,  relative  to  their  body  axis  (Nachtigall  et  al.,  1971; 

483  Esch  et  al.,  1975;  Taylor,  2001).  If  this  is  indeed  the  case,  then  the  only  way  to  increase  flight 

484  speed  -  and  to  continue  to  fly  horizontally  -  would  be  to  decrease  the  angle  of  elevation  of  the 

485  force  vector  by  pitching  the  thorax  downwards,  and  raising  the  abdomen  (Nachtigall  et  al., 
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486  1971).  Thus,  our  observed  changes  in  abdominal  pitch  may  reflect  a  part  of  this  underlying 

487  strategy  of  flight  control.  However,  exactly  how  changes  in  abdominal  elevation  can  affect 

488  changes  in  the  pitch  of  the  thorax  is  not  fully  understood,  although  a  model  examining  these 

489  changes  has  been  proposed  for  the  fruitfly  Drosophila  fZanker,  1988),  and  the  answer  would 

490  depend  partly  on  the  (as  yet  undetermined)  location  of  the  insect’s  centre  of  gravity  relative 

491  to  the  point  of  action  of  the  net  force  vector  generated  by  the  wings  (Taylor,  2001). 

492 

493  It  should  be  noted  that  the  responses  described  in  the  present  study  have  been  obtained  using 

494  an  ‘open  loop’  mode  of  experimentation  in  which  the  insect  is  responding  to  the  stimulus  that 

495  it  receives,  but  is  unable  to  influence  the  stimulus.  Our  experiments  reveal  that  body  posture 

496  is  strongly  influenced  by  the  visual  stimulus,  but  they  do  not  tell  us  whether  an  unrestrained 

497  bee  in  free  flight  controls  its  flight  speed  (and  therefore  the  image  speed  that  it  experiences) 

498  by  adjusting  its  body  posture.  Further  studies,  conducted  using  a  ‘closed  loop’  paradigm  in 

499  which  the  tethered  insect  is  able  to  freely  adjust  the  pitch  of  its  body,  and  where  the  measured 

500  thrust  and  lift  are  used  to  control  the  motion  of  the  visual  panorama  as  well  as  the  airflow 

501  appropriately,  should  provide  a  more  complete  understanding  of  the  phenomenon  that  we 

502  have  uncovered  here. 

503 

504 

505 

506 
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628  Figure  legends 

629 

630  Figure  1.  The  physical  and  virtual  reality  setup  for  tethered  hee  experiments.  (A)  Honeybees 

631  were  tethered  to  the  horizontal  end  of  an  L-shaped  metallic  rod.  Flight  was  initiated  hy 

632  lowering  a  supporting  platform.  (B)  Each  honeybee  was  placed  in  the  centre  of  the  arena 

633  composed  of  four  monitors.  A  video  camera  was  used  to  film  a  side  view  of  the  bee  in  flight. 

634  (C)  The  virtual  scene  consisted  of  an  infinite  tunnel  composed  two  walls  lined  with  vertically 

635  oriented  red  and  white  sinusoidal  gratings.  The  honeybee  was  at  a  virtual  distance  of  X  units 

636  from  each  of  the  tunnel  walls,  and  experienced  virtual  forward  motion  at  a  speed  of  Y  units 

637  per  display  frame. 

638 

639  Figure  2.  Illustration  of  the  measurement  of  the  streamlining  response,  defined  as  the 

640  orientation  of  the  abdomen  relative  to  the  thorax.  The  upper  row  of  sketches  illustration  the 

641  definition  of  the  response,  and  the  lower  row  show  two  examples  of  how  this  is  used  to 

642  measure  the  response.  Details  in  text. 

643 

644  Figure  3.  Comparison  of  streamlining  responses  induced  by  forward  image  motion  (FIM), 

645  no  image  motion  (NIM)  and  backward  image  motion  (BIM).  The  grey  staircase  displays  the 

646  progressively  ascending  stimulus  speed  (see  Methods).  Each  data  point  shows  mean  abdomen 

647  angle  ±  s.e.m.  Note  the  comparatively  small  variability  in  the  responses  when  bees  were 

648  experiencing  forward  image  motion. 

649 

650  Figure  4.  Streamlining  responses  of  tethered  bees  in  experiments  using  various 

651  configurations  of  two  monitors  (colored  symbols),  compared  with  the  four-monitor 

652  configuration  (black  symbols).  All  flights  were  presented  with  the  ascending  image  motion 

653  speed  protocol,  as  indicated  by  the  grey  staircase.  Each  data  point  shows  mean  abdomen 

654  angle  ±  s.d.  Data  comprised  only  from  bees  that  experienced  all  5  monitor  configurations 

655  (n=7). 

656 

657  Figure  5.  Control  experiments  showing  the  variation  of  the  streamlining  response  with  time 

658  in  a  protocol  in  which  the  stimulus  speed  was  held  constant  throughout  the  experiment  at  355 

659  °/s  or  at  526  °/sec.  There  is  no  systematic  variation  of  the  mean  response  with  time, 

660  indicating  that  response  fatigue  is  not  a  factor  influencing  the  results  in  the  experiments  of 

661  figure  4.  Each  data  point  shows  the  mean  abdomen  angle  ±  s.e.m. 
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662 

663  Figure  6.  Effect  of  stimulus  contrast  on  the  streamlining  response.  Variation  of  abdomen 

664  orientation  with  stimulus  speed  for  five  different  contrasts,  as  shown  in  the  index.  Each  data 

665  point  shows  the  mean  abdomen  angle  ±  s.e.m.  (n=6). 
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Appendix 


Table  Al.  Statistical  results  for  comparisons  made  between  the  various  two- 
monitor  configurations. 

_ Mauchly's  Test  of  Sphericity‘s _ 


Within  Subjects 

Epsilon® 

Effect 

Mauchly's 

W 

Approx.  Chi- 

Square 

df 

Sig. 

Greenhouse- 

Geisser 

Huynh- 

Feldt 

Lower- 

bound 

MonConfig 

.195 

7.713 

5 

.181 

.585 

.808 

.333 

OFR 

.000 

35 

.193 

.246 

.125 

MonConfig  *  OFR 

.000 

299 

.155 

.437 

.042 

Tests  the  null  hypothesis  that  the  error  covariance  matrix  of  the  orthonormalized  transformed  dependent 
variables  is  proportional  to  an  identity  matrix. 

a.  May  be  used  to  adjust  the  degrees  of  freedom  for  the  averaged  tests  of  significance.  Corrected  tests  are 
displayed  in  the  Tests  of  Within-Subjects  Effects  table. 

b.  Design:  Intercept 

Within  Subjects  Design:  MonConfig  +  OFR  +  MonConfig  *  OFR 


Tests  of  Within-Subjects  Effects 


Source 

Type  III  Sum 

of  Squares 

df 

Mean  Square 

LL 

MonConfig 

Sphericity  Assumed 

1029.372 

3 

343.124 

.437 

.729 

Greenhouse-Geisser 

1029.372 

1.756 

586.308 

.437 

.632 

Huynh-Feldt 

1029.372 

2.424 

424.593 

.437 

.691 

Lower-bound 

1029.372 

1.000 

1029.372 

.437 

.533 

Error(MonConfig) 

Sphericity  Assumed 

14138.202 

18 

785.456 

Greenhouse-Geisser 

14138.202 

10.534 

1342.135 

Huynh-Feldt 

14138.202 

14.546 

971.949 

Lower-bound 

14138.202 

6.000 

2356.367 

OFR 

Sphericity  Assumed 

42785.115 

8 

5348.139 

12.910 

.000 

Greenhouse-Geisser 

42785.115 

1.541 

27768.953 

12.910 

.003 

Huynh-Feldt 

42785.115 

1.970 

21716.948 

12.910 

.001 

Lower-bound 

42785.115 

1.000 

42785.115 

12.910 

.011 

Error(OFR) 

Sphericity  Assumed 

19885.242 

48 

414.276 

Greenhouse-Geisser 

19885.242 

9.245 

2151.030 

Huynh-Feldt 

19885.242 

11.821 

1682.231 

Lower-bound 

19885.242 

6.000 

3314.207 

MonConfig  *  OFR 

Sphericity  Assumed 

2954.672 

24 

123.111 

.811 

.718 

Greenhouse-Geisser 

2954.672 

3.714 

795.521 

.811 

.524 

Huynh-Feldt 

2954.672 

10.499 

281.428 

.811 

.624 

Lower-bound 

2954.672 

1.000 

2954.672 

.811 

.402 

Error(MonConfig*OFR)  Sphericity  Assumed 

21856.786 

144 

151.783 

Greenhouse-Geisser 

21856.786 

22.285 

980.793 

Huynh-Feldt 

21856.786 

62.993 

346.971 

Lower-bound 

21856.786 

6.000 

3642.798 

Table  A2.  Statistical  results  for  comparisons  made  between  the  two-  and  four- 
monitor  configurations. 


Mauchly's  Test  of  Sphericity*^ 


Mauchly's 

W 

Approx.  Chi- 

Square 

df 

Sig. 

Epsilon®  I 

Greenhouse- 

Geisser 

Huynh- 

Feldt 

Lower- 

bound 

MonConfig 

.173 

7.739 

9 

.589 

.625 

1.000 

.250 

OFR 

.000 

35 

.182 

.226 

.125 

MonConfig  *  OFR 

.000 

527 

.128 

.439 

.031 

Tests  the  null  hypothesis  that  the  error  covariance  matrix  of  the  orthonormalized  transformed  dependent  variables 
is  proportional  to  an  identity  matrix. 

a.  May  be  used  to  adjust  the  degrees  of  freedom  for  the  averaged  tests  of  significance.  Corrected  tests  are 
displayed  in  the  Tests  of  Within-Subjects  Effects  table. 

b.  Design:  Intercept 

Within  Subjects  Design:  MonConfig  +  OFR  +  MonConfig  *  OFR 


Tests  of  Within-Subjects  Effects 


Source 

Type  III  Sum 

of  Squares 

df 

Mean  Square 

F 

MonConfig 

Sphericity  Assumed 

59954.038 

4 

14988.509 

19.617 

.000 

Greenhouse-Geisser 

59954.038 

2.501 

23968.601 

19.617 

.000 

Huynh-Feldt 

59954.038 

4.000 

14988.509 

19.617 

.000 

Lower-bound 

59954.038 

1.000 

59954.038 

19.617 

.004 

Error(MonConfig) 

Sphericity  Assumed 

18336.969 

24 

764.040 

Greenhouse-Geisser 

18336.969 

15.008 

1221.801 

Huynh-Feldt 

18336.969 

24.000 

764.040 

Lower-bound 

18336.969 

6.000 

3056.162 

OFR 

Sphericity  Assumed 

45500.309 

8 

5687.539 

14.298 

.000 

Greenhouse-Geisser 

45500.309 

1.457 

31227.072 

14.298 

.003 

Huynh-Feldt 

45500.309 

1.805 

25209.218 

14.298 

.001 

Lower-bound 

45500.309 

1.000 

45500.309 

14.298 

.009 

Error(OFR) 

Sphericity  Assumed 

19093.636 

48 

397.784 

Greenhouse-Geisser  19093.636  8.742  2184.008 

Huynh-Feldt  19093.636  10.829  1763.122 

Lower-bound  19093.636  6.000  3182.273 


MonConfig  *  OFR 

Sphericity  Assumed 

6165.922 

32 

192.685 

1.379 

.098 

Greenhouse-Geisser 

6165.922 

4.100 

1503.859 

1.379 

.270 

Huynh-Feldt 

6165.922 

14.053 

438.751 

1.379 

.182 

Lower-bound 

6165.922 

1.000 

6165.922 

1.379 

.285 

Error(MonConfig*OFR)  Sphericity  Assumed  I  26836.069  192  139.771 


Greenhouse-Geisser  26836.069  24.600  1090.880 

Huynh-Feldt  26836.069  84.320  318.264 

Lower-bound _ 26836.069  6.000 _ 4472.678 

_ Pairwise  Comparisons _ 


(1)  MonConfig  (J)  MonConfig 

Mean 

Std.  Error 

Sig." 

95%  Confidence  Interval  for 

Difference® 

Difference  (l-J) 

Lower  Bound 

Upper  Bound 

1  2 

2.648 

6.919 

1.000 

-27.218 

32.515 

3 

-2.548 

2.243 

1.000 

-12.228 

7.133 

4 

1.979 

4.276 

1.000 

-16.478 

20.436 

5 

-33.673* 

4.424 

.003 

-52.771 

-14.575 

2  1 

-2.648 

6.919 

1.000 

-32.515 

27.218 

3 

-5.196 

5.617 

1.000 

-29.442 

19.050 

4 

-.670 

5.595 

1.000 

-24.821 

23.482 

5 

-36.321* 

5.559 

.006 

-60.317 

-12.326 

3  1 

2.548 

2.243 

1.000 

-7.133 

12.228 

2 

5.196 

5.617 

1.000 

-19.050 

29.442 

4 

4.526 

3.948 

1.000 

-12.515 

21.568 

5 

-31.125* 

3.936 

.002 

-48.117 

-14.133 

4  1 

-1.979 

4.276 

1.000 

-20.436 

16.478 

2 

.670 

5.595 

1.000 

-23.482 

24.821 

3 

-4.526 

3.948 

1.000 

-21.568 

12.515 

5 

-35.652* 

5.194 

.005 

-58.073 

-13.230 

5  1 

33.673* 

4.424 

.003 

14.575 

52.771 

2 

36.321* 

5.559 

.006 

12.326 

60.317 

3 

31.125* 

3.936 

.002 

14.133 

48.117 

4 

35.652* 

5.194 

.005 

13.230 

58.073 

Based  on  estimated  marginal  means 
a.  Adjustment  for  multiple  comparisons:  Bonferroni. 
*.  The  mean  difference  is  significant  at  the  .05  level. 


Table  A3.  Statistical  results  for  comparisons  made  between  the  four  non-zero 
contrast  conditions. 


Mauchly's  Test  of  Sphericity‘s 


Mauchly's 

W 

Approx.  Chi- 

Square 

df 

Sig. 

Epsilon® 

Greenhouse- 

Geisser 

Huynh- 

Feldt 

Lower- 

bound 

Contrast 

.608 

1.853 

5 

.873 

.745 

1.000 

.333 

OFR 

.000 

35 

.275 

.499 

.125 

Contrast  *  OFR 

.000 

299 

.104 

.214 

.042 

Tests  the  null  hypothesis  that  the  error  covariance  matrix  of  the  orthonormalized  transformed  dependent  variables 
is  proportional  to  an  identity  matrix. 

a.  May  be  used  to  adjust  the  degrees  of  freedom  for  the  averaged  tests  of  significance.  Corrected  tests  are 
displayed  in  the  Tests  of  Within-Subjects  Effects  table. 

b.  Design:  Intercept 

Within  Subjects  Design:  MonConfig  +  OFR  +  MonConfig  *  OFR 


Tests  of  Within-Subjects  Effects 


Source 

Type  III  Sum 

of  Squares 

df 

Mean  Square 

F 

Sig. 

Contrast 

Sphericity  Assumed 

4084.904 

3 

1361.635 

1.467 

.263 

Greenhouse-Geisser 

4084.904 

2.234 

1828.600 

1.467 

.274 

Huynh-Feldt 

4084.904 

3.000 

1361.635 

1.467 

.263 

Lower- bound 

4084.904 

1.000 

4084.904 

1.467 

.280 

Error(Contrast) 

Sphericity  Assumed 

13921.102 

15 

928.073 

Greenhouse-Geisser 

13921.102 

11.169 

1246.351 

Huynh-Feldt 

13921.102 

15.000 

928.073 

Lower- bound 

13921.102 

5.000 

2784.220 

OFR 

Sphericity  Assumed 

40476.934 

8 

5059.617 

13.218 

.000 

Greenhouse-Geisser 

40476.934 

2.197 

18423.584 

13.218 

.001 

Huynh-Feldt 

40476.934 

3.989 

10146.225 

13.218 

.000 

Lower- bound 

40476.934 

1.000 

40476.934 

13.218 

.015 

Error(OFR) 

Sphericity  Assumed 

15310.864 

40 

382.772 

Greenhouse-Geisser 

15310.864 

10.985 

1393.786 

Huynh-Feldt 

15310.864 

19.947 

767.585 

Lower- bound 

15310.864 

5.000 

3062.173 

Contrast  *  OFR 

Sphericity  Assumed 

1688.328 

24 

70.347 

.599 

.927 

Greenhouse-Geisser 

1688.328 

2.486 

679.146 

.599 

.599 

Huynh-Feldt 

1688.328 

5.137 

328.632 

.599 

.705 

Lower- bound 

1688.328 

1.000 

1688.328 

.599 

.474 

Error(Contrast*OFR)  Sphericity  Assumed  14097.780  120  117.482 

Greenhouse-Geisser  14097.780  12.430  1134.193 

Huynh-Feldt  14097.780  25.687  548.825 

Lower-bound  14097.780  5.000  2819.556 


Table  A4.  Statistical  results  for  comparisons  made  between  control  (zero)  and 
four  contrast  conditions. 


Mauchly's  Test  of  Sphericity‘s 


Mauchly's 

W 

Approx.  Chi- 

Square 

df 

Sig. 

Epsilon® 

Greenhouse- 

Geisser 

Huynh- 

Feldt 

Lower- 

bound 

Contrast 

.266 

4.527 

9 

.891 

.684 

1.000 

.250 

OFR 

.000 

35 

.294 

.572 

.125 

Contrast  *  OFR 

.000 

527 

.077 

.156 

.031 

Tests  the  null  hypothesis  that  the  error  covariance  matrix  of  the  orthonormalized  transformed  dependent  variables 
is  proportional  to  an  identity  matrix. 

a.  May  be  used  to  adjust  the  degrees  of  freedom  for  the  averaged  tests  of  significance.  Corrected  tests  are 
displayed  in  the  Tests  of  Within-Subjects  Effects  table. 

b.  Design:  Intercept 

Within  Subjects  Design:  MonConfig  +  OFR  +  MonConfig  *  OFR 


Tests  of  Within-Subjects  Effects 


Source 

Type  III  Sum 

of  Squares 

df 

Mean  Square 

F 

Sig. 

Contrast 

Sphericity  Assumed 

235379.647 

4 

58844.912 

53.239 

.000 

Greenhouse-Geisser 

235379.647 

2.737 

85991.499 

53.239 

.000 

Huynh-Feldt 

235379.647 

4.000 

58844.912 

53.239 

.000 

Lower-bound 

235379.647 

1.000 

235379.647 

53.239 

.001 

Error(Contrast) 

Sphericity  Assumed 

22105.838 

20 

1105.292 

Greenhouse-Geisser 

22105.838 

13.686 

1615.190 

Huynh-Feldt 

22105.838 

20.000 

1105.292 

Lower-bound 

22105.838 

5.000 

4421.168 

OFR 

Sphericity  Assumed 

28356.376 

8 

3544.547 

10.672 

.000 

Greenhouse-Geisser 

28356.376 

2.352 

12053.893 

10.672 

.002 

Huynh-Feldt 

28356.376 

4.576 

6196.923 

10.672 

.000 

Lower-bound 

28356.376 

1.000 

28356.376 

10.672 

.022 

Error(OFR) 

Sphericity  Assumed 

13285.319 

40 

332.133 

Greenhouse-Geisser 

13285.319 

11.762 

1129.480 

_ Huynh-Feldt 

13285.319 

22.879 

580.667 

Lower-bound  13285.319  5.000  2657.064 


Contrast  *  OFR  Sphericity  Assumed 

15339.390 

32 

479.356 

4.126 

.000 

Greenhouse-Geisser 

15339.390 

2.450 

6261.070 

4.126 

.037 

Huynh-Feldt 

15339.390 

4.980 

3080.054 

4.126 

.007 

Lower-bound 

15339.390 

1.000 

15339.390 

4.126 

.098 

Error(Contrast*OFR)  Sphericity  Assumed  I  18590.789  160  116.192 


Greenhouse-Geisser  18590.789  12.250  1517.638 

Huynh-Feldt  18590.789  24.901  746.583 

Lower-bound _  18590.789 1  5.000  |  3718.158 

_ Pairwise  Comparisons _ 


(1)  Contrast  (J)  Contrast 

Mean 

Std.  Error 

Sig." 

95%  Confidence  Interval  for 

Difference® 

Difference  (l-J) 

Lower  Bound 

Upper  Bound 

1  2 

-71.907* 

8.769 

.004 

-113.766 

-30.048 

(C=0)  3 

-66.671* 

6.882 

.002 

-99.519 

-33.822 

4 

-76.146* 

6.158 

.001 

-105.540 

-46.752 

5 

-77.961* 

6.399 

.001 

-108.507 

-47.415 

2  1 

71.907* 

8.769 

.004 

30.048 

113.766 

(C=0.29)  3 

5.237 

7.055 

1.000 

-28.441 

38.915 

4 

-4.239 

4.111 

1.000 

-23.862 

15.384 

5 

-6.054 

6.348 

1.000 

-36.356 

24.249 

3 

'1 

66.671* 

6.882 

.002 

33.822 

99.519 

(C=0.36)  ^ 

-5.237 

7.055 

1.000 

-38.915 

28.441 

4 

-9.475 

6.308 

1.000 

-39.585 

20.634 

5 

-11.290 

5.040 

.752 

-35.349 

12.769 

4 

'1 

76.146* 

6.158 

.001 

46.752 

105.540 

(C=0.54)  ^ 

4.239 

4.111 

1.000 

-15.384 

23.862 

3 

9.475 

6.308 

1.000 

-20.634 

39.585 

5 

-1.815 

5.837 

1.000 

-29.675 

26.045 

5 

'1 

77.961* 

6.399 

.001 

47.415 

108.507 

(C=0.72)  ^ 

6.054 

6.348 

1.000 

-24.249 

36.356 

3 

11.290 

5.040 

.752 

-12.769 

35.349 

4 

1.815 

5.837 

1.000 

-26.045 

29.675 

Based  on  estimated  marginal  means 
a.  Adjustment  for  multiple  comparisons:  Bonferroni. 
*.  The  mean  difference  is  significant  at  the  .05  level. 


